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1,   Introduction 

iMo   theories  of  the  propagation  of  spherical  shocks,  having 
a  similar  approach,  have  been  put  forward  Independently.   Cne 
theory  was  developed  by  Klrkwood  and  Bethe  in  1942,  but  then  it 
was  contained  in  a  confidential  report,  and  the  first  published 
accovmt  appeared  (without  full  details)  in  1948  in  R»  H.  Cole's 
book  "Underwater  Exploslom" ,  Unaware  of  this  work,  the  author 
obtaimed  a  theory,  early  In  1951  #  which  a  recent  iBuvestlgatlon 
has  shown  to  have  much  in  common  with  the  Klrkwood-Bethe  theory. 
Klrkwood  and  Bethe  considered  the  problem  of  an  explosion  in 
water  and  found  how  the  head  shock  propagates;  it  is  relevant  to 
note  here  that  their  results  are  intended  to  hold  tme    (for  water) 
even  if  the  shock  is  not  weak.   On  the  other  hand,  the  author 
found  a  theory  for  the  propagation  of  weak  spherical  shocks  (in 
any  medium)  but  obtained  a  description  of  the  whole  flow;  for 
example,  later  shocks  which  occur  in  addition  to  the  head  shock 
were  predicted  and  determined.   Thus  the  theories  overlap  in 
that  they  both  apply  to  a  weak  head  shock;  the  Klrkwood-Bethe 
theory  treats,  in  addition,  a  non-weak  head  shock  in  water  (an 
addition  of  much  practical  value),  whilst  the  author's  theory 
determines,  in  addition,  the  later  flows  and  additional  shocks. 

However,  Klrkwood  and  Bethe  made  various  assumptions,  the 
validity  of  which  was  not  clear,   l^ow,  when  their  theory  is  re- 
considered in  the  light  of  the  later  theory  of  weak  shocks  and 
derived  in  a  similar  way,  it  becomes  apparent  how  the  extension 
to  a  non-weak  head  shock  in  water  may  be  achieved  and  the  assiimp- 
tlona  involved  (although  they  still  cannot  be  completely  justi- 
fied) become  clear  and  appear  reasonable.   This  relation  between 
the  two  theories  is  brought  out  in  the  present  report.   First, 
in  Part  I  the  author's  theory  of  weak  shocks  is  described,  then 
in  Part  II  the  Klrkwood-Bethe  theory  is  presented  (but  with  the 
simpler  derivation  following  Part  I)  as  a  method  of  improving 
the  results  to  apply  to  a  non-weak  head  shock,  in  the  case  when 
the  medium  is  water.  Although  it  is  recom?nended  that  the  results 
now  be  vieved  in  this  way,  it  should  be  remembered  that  it  is 
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the  reverse  of  the  chronological  order  and  that  most  of  the  re- 
sults for  the  head  shock  in  part  I  were  also  obtained  (or  could 
have  been  obtained)  In  the  vjork  of  'rlrkwood-Bethe. 

Part  I.   The  Propagation  of  lueak  Shocks 
2.   Description  of  the  method 

The  propagation  of  weak  spherical  waves  and  shocks  will  be 
described  following  the  method  of  Fef ,  1,  in  which  the  shocks 
occurring  in  steady  supersonic  flow  past  an  axlsymrietrlcal  body 
were  determined.   The  direct  analogy  to  the  latter  problem  is 
the  propagation  of  cylindrical  shocks;  but  the  case  of  spherical 
symmetry  is  similar.   The  starting  point  of  this  method  is  the 
linearised  theory,  which  in  this  case  is  the  simple  acoustic 
theory.   The  linear  results,  although  extremely  valuable  for 
certain  purposes,  become  inadequate  to  describe  the  propagation 
after  some  time  since  they  do  not  give  the  correct  "spreading  out" 
of  the  wave;  the  linearized  theory  does  not  allow  for  the  varia- 
tion in  the  local  sovind  speed  or  for  the  connection  of  sound 
with  the  fluid,  and  as  time  goes  on  these  small  variations  ac- 
cumulate.  Moreover,  shocks,  which  play  an  important  part  in  the 
flow,  cannot  be  included  in  such  a  theory.   However,  these  defects 
can  be  corrected  to  give  a  valid  description  of  the  non-linear 
propagation;  this  is  the  basis  of  the  present  theory. 

A  hypothesis  is  made  that  the  Improvement  of  linearized  theory 
may  be  achieved  simply  by  replacing  the  linearized  characteristics 
occurring  therein  by  a  more  (and  sufficiently)  accurate  approxi- 
mation to  the  exact  characteristics;  that  is,  it  is  claimed  that 
the  linear  results  show  the  correct  variation  of  the  quantities 
along  the  characteristics  but  do  not  have  the  correct  curves 
(In  the  (r,t)  plane)  for  the  characteristics.   Since  character- 
istics are  the  paths  of  wavelets,  this  hypothesis  means  that  the 
propagation  speed  of  linearized  theory,  which  is  c   (the  soiind 
speed  in  the  undisturbed  fluid),  is  not  a  satisfactory  approxi- 
mation to  the  exact  propagation  speed  c+u  (where  c  Is  the  local 
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sound  speed  and  u  is  the  fluid  velocity)  and  must  be  corrected. 
Hence  the  first  failing  of  the  linear  results  is  immediately 
rectified.   The  other,  concerning  shocks,  also  disappears  because 
now  the  characteristics  are  not  fixed  parallel  curves;  in  cer- 
tain regions  they  run  tof ether  in  compression,  predicting  shocks, 
and  the  requisiteshocks  (to  prevent  the  characteristics  from 
intersecting)  are  determined  from  the  Pankine-Hugoniot  shock 
conditions.   In  connection  with  shocks  it  is  necessary  to  point 
out  that  the  linearized  results  and  the  theory  which  is  deduced 
from  them  neglect  chajiges  In  the  entropy  (so  that  the  pressure 
is  a  function  only  of  the  density).  This  is  valid,  however, 
because  the  entropy  Jximp  at  a  shock  is  of  third  order  in  the 
strength  of  the  shock;  hence  for  the  weak  shocks  considered  here, 
the  entropy  changes  may  be  neglected  in  obtaining  the  first  order 
results  for  the  flow  behind,   Ftirther  comments  on  this  point  may 
be  found  in  Ref.  1,  p.  303. 

In  addition  to  the  above  physical  considerations,  the  hypo- 
thesis is  substantiated  in  many  ways  (they  are  discussed  in 
Ref,  1  and  are  not  repeated  here)  and  its  validity  should  not 
be  considered  to  be  in  doubt. 

It  is  supposed  for  definiteness  that  some  initial  particle 
path  r  =  a(t),  say  is  given  (where  r  and  t  are  the  distance 
from  the  center  of  syninetry  and  t  is  the  time),  and  a  description 
of  the  propagation  outside  this  path  is  required.   The  problem 
could  be  set  up  in  other  ways  but  this  one  is  found  to  be  very 
convenient.   For  example,  in  an  underwater  explosion  problem, 
a(t)  would  be  the  radius  of  the  ^.as  sphere.   Alternatively, 
r  =  a(t)   could  be  taken  as  a  prescribed  variation  v/ith  time  of 
the  radius  of  a  "spherical  piston."  Since  the  theory  only 
applies  to  weak  disturbances  it  is  assumed  that  the  velocity 
a'(t)  at  the  particle  path  is  small  compared  to  the  undisturbed 
sound  speed  c  ;  the  small  quantity  5  which  measures  the  smallness 
of  the  disturbance  is  defined  as  max^a'(t)/c  f  .   If  the  explo- 
sion is  not  weak  the  present  theory  does  not  apply  in  the  initial 
stages.   But,  a  real  explosion  always  decays  so  that  the  dis- 
turbance will  be  weak  at  a  sufficient  distance  from  the  center. 
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In  this  region  a  particle  path  could  be  taken  (the  data  for  It 
being  obtained  by  numerical  methods,  for  example)  and  the  theory- 
applied  to  find  the  subsequent  propagation  to  infinity. 

3 •   Llnearl zed  theo ry  and  its  Impro vement 

The  linearized  (acoustic)  theory  of  spherical  waves  is  ob- 
tained by  neglecting  second  and  higher  order  terms  in  the  equa- 
tions of  motion;  the  following  results  are  well-known.   If  a 
potential  cj  is  introduced  such  that  u  =  b^/br,    then  4  satisfies 


(1)     ipp  -  I  "tr  =  -T  Kt      - 

and  the  pressure  p  Is  given  in  terms  of  4  by 

(2)       P  -  Po  =  -Po  '^t    ' 

where  p  and  p  are  the  undisturbed  values  of  p  and  the  density 
o      o 

p.   The  equation  of  state  is  assumed  to  be 

(3)  p  =  KiP^  -  ko     * 

where  y,  }r  ,  |(,  are  constant'';  usually  gases  are  considered  and 
k     Is  zero,  but  for  water  "^    is  Important  (its  value  being  ap- 
proximately 3,000  atmospheres).   Then  the  local  sound  speed  c  Is 
found  by  approrlmatlnf;  the  Bernculll  equation 

c^   ^  1  j2  ^  .   _  °o 

to  give 

(i^)       c  ^  c^  -  ^  (Y-l)4tA^    . 

The  solution  of  (1)  representing  sn  outgoing  wave  Is 

4  =  -f(t-r/c^)/r   ; 


In  general  k     and  ^  are  functions  of  the  entropy,  but  entropy 
changes  are  negligible  here. 
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c  =  c 
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P  -  Pg  =  P0 
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These  results  are  corrected,  as  explained  In  q2,    by  replacing 
the  linearized  characteristic  variable   t-r/c  ,  which  appears 
in  them  by  c'(ts,r)  where  the  curves   7r  =  constant  are  exact 
characteristic  curves  (or  at  least  sufficiently  good  approxima- 
tions to  them) .   The  definition  of  T    is  not  complete  since  it 
is  only  required  to  be  constant  on  a  characteristic;  it  is  made 
unique  by  defining  it  to  be  the  value  of  t  where  the  character- 
istic meets  the  particle  path  r  =  a(t).   Kodlfying  the  linear 
theory  in  this  way,  and  at  the  same  time  replacing  the  function 
f  by   f   P('Ci)dt',  the  results  become 


(6) 
(7) 
(8) 


I\ow  ^(tjr)  Is  determined  from  the  condition  that,  on  a  charac- 
teristic T=   constant,  dr/dt  =  u+c ;  therefore,  from  (6)  and  (7), 
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vjhere     k  =  ^{y+l)/<i^»      Equation   (8)    is   intef;r3ted,   under  the 
condition  that     t  =T     when     r  =  a{r),    to  give 

(10)     t  =  £^i  -  £KP(r)log  5^  -  i2— 2 /_^.i).r 

o 

Of  course,  since  the  expressions  (6)  and  (?)  are  not  exact  and 
in  (9)  terir.s  of  second  and  higiier  order  in  u  and  c-c  are  neg- 
lected, c'(t,r)  is  still  not  the  exact  characteristic  variable; 
however,  in  contrast  to  the  characteristics  of  linearized  theory, 
(IC)  is  a  valid  approximation  to  the  exact  characteristic  carves, 
because  the  tern.s  neglected  in  it  really  are  small  compared  to 
those  retained,  Llneariiied  theory,  for  example,  neglected  the 
term  -2kP('2:)log  af"\      ^^^  comparison  with  t-  -^  on  the  basis 

that  ?  is  small;  but  this  term  becomes  infinitely  large  com- 
pared to   t-  —  both  when  r  — >  oo  and  vjhen  t-  —  is  near 

c  o 

zero  I   Equations  (6),  (7),  (8)  and  (10)  form  the  corrected  solu- 
tion of  the  flowa 

The  arbitrary  function  F{Z)   which  occurs  in  the  solution 
must  be  detennined  from  the  given  conditions  on  the  initial 
curve   r  =  a(t)o   Since  r  =  a(t)   is  a  particle  path,  u  =  a'(t) 
on  it;  hence,  from  (6)  remembering  that  on  r=a(t),t=?. 


(^1)       ^'t^)  =^4m^-^ 


/  F(t')dt» 


a"(t) 

An  integrating  factor  for  this  equation  is  c  a^x,  where 


tj.(t)  =  exp  I  c^dt/a(t),  whence 

J         o 

(12)       /  F(t')dt«  =  -t|y  ,/"  a(7)a'(7)^(v^)dV 


2 


ft 
(13)      F(t)  =  c^a(t)a«(t)  -  3^(^)°(^)  I  a(7)a•(7)^i(7)dV 

If,  in  addition  to  knowing  the  particle  path,  the  value  of  p  on 
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It  Is  also  known  as  a  function  of  time,  so  that  P-P-  =  ?(t) 
say,  then  P  Is  detemlned  more  simply  from  (8)  as 

(li^)      P(t)  =  p;^a(t)P{t) 

In  the  special  case  when  a(t)  =  0  for  t  <  0,  and  a"(t)  Is 
continuous  for  t  >  0,  a  simpler  expression  than  (13)  for  F  may 
be  obtained!  the  second  term  in  (11)  dominates  the  first  (since 
the  ratio  of  the  first  term  to  the  second  is 
c(a(t)/(c  t))  =  0(a«(t)/c  )  =»  0(6))  and  P  may  be  taken  as 

{1$)  P(t)  =  ^  (a^aM   . 

Por  a  real  explosion,  p  returns  to  p  and  u  returns  to  C 
after  an  infinite  timej  hence  from  (6)  and  (8)  it  is  clear  that 
i^(f)  _>  0  as  T->  0     and  that 
p  oo 

(16)  /     p(r)dr  =  0 

Jo 
Later,  this  result  will  be  shown  to  have  very  Important  conse- 
quences.  It  may  also  be  deduced  from  (13)  by  assuming  that  a(t) 
is  bovmded  and  that  a»(t)  ->  C  as  t  —>  oo   (assumptions  which 
would  be  true  for  a  real  explosion)  but  the  simple  arg\iment  above 
is  more  instructivee   In  general,  the  fluid  is  pushed  outward 
suddenly  so  that  a'(t)  is  discontinuous  and  positive  at  t  =  0; 
hence  from  (13)  P(T)  is  discontinuous  and  positive  at    =  0. 
The  simplest  function  F(t)   would  be  when  a'(t)  has  an  Initial 
positive  value  c  e,  say  (where  e  is  small  and  0(6)),  and  tends 

monotonlcally  to  zero  as  t  — >  oo  .   Then  F{T)   has  the  initial 

2 

positive  value  c  a  e,  where  a  =  a(0),  and  is  of  the  form  shown 

o  o  '        o 

in  figure  1,  the  aa:'eas  are  as  above  and  below  the  axis  being 
equal  by  (16).   '^  P 


r- ^-l"-      - -,^^-..-.r 


1 


.-rr^" 


Figure  1. 
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l^.o      The  shocks 

The  solution  obtained  In  the  last  section  is  not  yet  com- 
plete.  The  characteristics  now  depend  on  the  solution;  where 
P'(t^)  <   0,  the  characteristics'"'  ( IC )  diverge  in  expansion,  and 
where   F'(T)  >  0   they  run  together  and  eventually  overlap. 
The  extreme  cases  of  "centered"  expansion  or  compression,  with 
a  fan  of  characteristics  from  a  single  point  on  r  =  a(t),  occur 
when  ?{T),   due  to  a  discontinuity  in  a»(t),  is  discontinuous 
(since  in  that  case  a  whole  range  of  values  of  P,  and  hence  a 
fan  of  characteristics  (10),  is  found  for  the  sarae  value  of  r )  • 
If  the  characteristics  were  allowed  to  overlap,  the  solution 
would  cease  to  be  single-valued  and  there  would  be  a  fold  in 
the  flow  plane;  a  shock  must  be  fitted  in  such  regions  so  that 
the  ch'iracteristlcs  meet  the  shock  and  are  cut  off  by  it  before 
they  meet  each  other.   Thus  shock  formation  is  predicted  in  the 
regions  where   ?'(f)  >  0  or  P(Tr)  has  a  positive  discontinuity; 
In  the  latter  case  the  shock  is  formed  Immediately  at  the  part- 
icle path  r  =  a(t)   (as  must  be  the  case  since  a'(t)  would  be 
discontinuous  there).   For  the  problem  of  interest  in  explosions, 
F((r)  has  been  seen  to  have  a  positive  discontinuity  at;   C  =  0; 
hence  there  will  be  a  shock.   This  head  shock  is  now  determined. 

( 1 )   The  head  shock 

The  Fanklne-Rugonlot  shock  conditions,  expressing  the  laws 
of  conservation  of  mass,  momentum,  and  energy  through  the  shock, 
may  be  written  (assiiming  that  the  equation  of  state  is  (3)  with 
U-t   ki  functions  of  the  entropy, 

(17)      ^  =  7^C   , 


For  many  purposes  the  term  c"*  /  F('r  •  )dt»  (a~  -r"  )   may  be 


o 

o 


neglected  since  it  is  small.   Certainly  questions  of  converg- 
ence or  divergence  of  characteristics  etc.  may  be  decided 

from  the  term  2kF('f)log  F  y      alone. 


.'f  :* 
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■r       ':   :   '  . 
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;  V'V  : 


-     i  ■■: . 


(16)     n2  =  c2  (i  .  |±1  . 

(19) 

where  U  is  the  shock  velocity  and  £  Is  defined  by 
(2C,     «  =  --5^   • 

In  the  case  uhen  ^q  =  0     In  (3)»  f  is  the  strenrth  (p-p  )/p  , 
otherwise  It  Is  the  strength  raultiplled  by  Pq/CPo'^'Ko)  •  '^°^'   *^® 
first  order  theory,  second  and  hlrher  order  temc  in  i   m-^y  be 
neglected,  so 

(21)  u  =  !2i   ,   u  =  c„  (l  *  Jgi  c)   .   c  =  o„(l  ^  1^  e)   , 

(At  the  sane  time  it  may  be  renirl-'ed  that  the  Jump  in  entrony 
at  the  shock  is  small  and  C(£,  );  this  is  the  reason  that  the 
change  of  entrop:/  J^^siy  be  neglected  In  deriving  tae  solution  in 
Section  3»)  Clearly  the  shoclc  conditions  furnish  tuo  relations: 
one  relation  beti;een  the  flovj  quantities  behind  the  shock,  and 
one  relation  to  determine  U.   It  is  convenient  to  choose  the.-n, 
from  (21),  as 

(22)  ^  =  :pr  ^°-%^  *  '^t^^J    » 

and 

(23)  U  =  I  (u+c+c^)  +  0(^^) 

(Actually  the  error  in  (22)  is  0(<.^).) 

Equation  (22)  leads  t"  an  important  point  vjhich  also  arises 
in  flndinr  the  shock  from  (23) «   lo  satisfy  {22),    it  is  necessary 
from  (6)  and  (7)  that,  at  the  shock,  the  second  term  in  (o)  is 
of  smaller  order  than  the  first;  that  is,  since  the  ratio  of  the 
terms  is  0(c  t/r),  it  is  required  that  c  r/r  is  small  on  the 
shock.   This  is  so.   I'or,  on  the  shock,  a  characteristic  of  the 


*-   ■*    .'  -  . 
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*■      I 


r.iK  ^^ 
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disturbed  flow  meets  a  characteristic  of  the  undisturbed  flow 
ahead  (in  fact,  (23)  shows  that  in  the  (r,t)  plane  the  shock 
bisects  the  angle  between  them);  on  the  undisturbed  character- 
istic, t-(r-a  )/c   is  equal  to  a  negative  constant,  therefore 


in  (IC) 


(2M   -^-^  *  21cP(r)log  ^  .  i^_2 (a^  -  ?)  >  ^ 


o 

at  the  shock.  i;ow  (a(r)-a^)/c^  and  ^"^  [   F(r«  )dr' {a""'"-r"*-'-) 

2  ° 
are  both  0(2:6);  thus,  since  P(f)  =  0(c^a  6),  (2tt)  can  only  be 

satisfied  if  r/a(r)  >  exp(i^c  T/a  S)')  where  K  is  a  certain 
constant.   -ence  c  2'/2*  Is  certainly  C(6),  and  is  actually  much 
smaller.   It  Is  proved,  then,  that  the  second  term  in  (6)  may 
be  neglected  at  the  shock  (in  agreement  with  (22)),  and,  more- 
over, that  its  contribution  to  (IC)  may  also  be  neglected. 
However,  avjay  from  the  snoci  its  presence  Is  extremely  important; 
for  example,  it  is  entirely  responsible  for  (lo)  which,  ac  will 
be  seen  later,  shews  that  at  least  one  other  shock  must  occur 
in  addition  to  the  head  shock. 
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Figure  2c 


If  the  equation  of  the  shock  is  assumed  to  be  t  =  (r-aj/c  -g(r) 


(7)  (with  the  second  term  in  (6)  neflected)  (23)  gives 


then  U  =  c  +c  g'(r)   to  a  fi_-st  approximation,  and  using  (6)  and 
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(25)  g'(r)  =  kF(r)r*"^   . 

Prom  the  argument  follovdng  (2l|.),  (IC)  may  be  approximated  to 

(26)  t  =  (r-a^)/cQ  -  2kP(r)log(r/a(r ))  +r   , 
at  the  shock;  hence  from  the  shock  equation 

(27)  g(r)  =  2kP(r)log(r/a(r))  -  T      . 

In  principle  equations  (25)  and  (27)  determine  the  f\inctlon  g(r) 
since  Z    could  be  eliminated;  it  is  more  convenient,  however,  to 
find  r  as  a  function  of  Z    i.e.  the  distance  at  which  the  shock 
cuts  the  characteristic  Z~   constant.   Llfferentiatin^  (27)  with 
respect  to  f  and  substltutinr  from  (26),  the  equation  for  v(7r) 
is  obtained: 

(28)  ^EFin  dr  ,  ,^p,^^)i,,  rr^  -  1  ^  ^^}X^^       , 


Sz   -"  ^«^'U';^oe,  ^^^  =  X  -^  ^j^ 


or, 


,z,,    ^p(.,xo,^j  =  .<„  .  !^Lm-(li 


fince   ;)t(  a«(r)  =  0/—   >'^M  =  0(5*^),  the  second  term  on  the 

^^^         V°o  °o  / 
right  hand  side  may  be  neglected;  hence 

(30)    kp2(r)iog  ^  =  r  F(rM<ir« 


It  may  be  pointed  out  that  in  (3C)  (and  in  (26)  in  the  neigh- 
borhood of  the  chock),  a(f)  can  be  approximated  by  a  in  the 
cases  when  a  ^0;  a(2r)  Itself  is  only  required  when  a  — >  0» 
(This  approximation  is  valid  because 


log  r/a^ ^   c^r/(a  6) ^  "  ^(S  ) 


from  the  discussion  following  (2I4.).   Hence  lo;?  r/a(2')  may  be 
approximated  by  log  r/a  .)  The  values  of  g  and  g'  are  found  by 
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substitution  in  (25)  and  (27).   The  strength  of  the  shock  could 
be  found  from  the  shock  conditions  but  is  most  easily  obtained 
from  (8), 

First  consider  the  initial  velocity  and  strength  of  the 

shock.   In  the  case  a  7^  C,  the  shock  starts  at  t  =  C,  r  =  a  , 

O  p  o 

with  initial  pressure  jump  p-p^  =  p^F(G)/a^  =  P^c^e  since 

p  "   '^      '^O     '^00     "^0  0 

?(C)  =  c  a  e,  where  ec  is  the  initial  velocity  of  the "piston 

00  o  ci") 

curve?';  the  initial  velocity  is  U^  =  c\l*  Mr-^l)e]    .   Both 

results  agree  with  the  results  for  plane  waves  (as  they  must 

since  the  spherical  nature  of  the  flow  ir  not  apcarent  initially). 

Ii.  the  case  a  =  C,  (1:?)  shows  that  P(?)  ~  2c-'^e-^r  for  small 

r;  hence  from  (3C) 

(31)  r  ^  c^er  exp  ^  (Y+l)'*'^e"^  ^    , 
from  (6),  the  initial  pressure  jump  is 

(32)  2p^c^e^eyp  j  ^(y+D'^e"^] 
and  from  (21) 

<33)      U^  =  °o  -^  I  (Y+l)Coe^exp  [-ir^ir^t'^] 

In  the  special  case  of  sphere  expanding  i-jith  uniform  velocity 
c  e,  ^^iZ)   =  2c^e^r   and  (31)-(33)  are  true  for  all  r;  there  is 
just  a  single  shock  of  constant  stren,  th  and  velocity  given  by 
(32)  and  (33)  •   In  cor.slderinf  the  uniformly  expanding  sphere 
Lifhthill  (i-vef.  2)  obtained  the  result  lor(u/cQ  -1)—  -(Y+l)''''"e'^ 
but  did  not  find  the  coefficient  of  the  exponential  in  (33). 

For  large  r,  (3C)  shors  that  ?'  is  near  to  T  ,    the  first 
zero  of  ?('^)  (a  zero  always  exists  by  (16)  for  a  real  explosion; 
for  an  artificial  example  such  as  the  uniformly  expanding  sphere 
considered  above  r  — >  00  as   r  — >  co  ) ,  and  (3C)  may  be  approxi- 
mated as  2kF(r)  =  b  log"-'-/^r/a(rQ),  where 

.2  _  Y+l  P° 


(31+)      h'^  =  ^J        F(r')d2'' 
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Therefore  g(r)  =  b  log  ^  (r/a(^  ))  -TqJ  that  is,  the  shock 


o 
arrives  at  the  distance  r  at  time 


(35)     b  Jioe  ^ 


earlier  than  the  (asymptotically)  undisturbed  characteristic 
t  =  (r-a^)/c 
Is  piven  by 


t  =  (r-a  )/c  +Z        (on   which  F  =  0).   From  (8),  the  pressure 


(36)      4  =  _j_o^  „  £X^    -j^__ 
Y   P„c^    ^    r  log  ' 


aTF"!" 


and  from  (21)  the  velocity  is 

c  b 
o 

2r  log-'-/^ 


(37)      U-  c^  + 


^r^ 


(11)   Later  shocks 

Condition  (16)  on  P  shows  that  there  must  always  be  a  sec- 
ond compressive  replon  (for  a  real  explosion)  In  which  ^'(r)  >  0 
and  the  characteristics  run  together  to  form  a  shock.   The  sim- 
plest case  will  be  an  F-curve  as  in  figure  1|  this  is  considered 
first.   Since  P('^)  Is  continuous,  the  distance  at  which  the 
characteristics  overlap  to  form  an  envelope  is  large.   For,  If 
(10)  is  written  as   t  =  h(2',r)  +2"   then  the  envelope  is 
^h/dr  +  1  =  0;  all  the  terms  in  'bh/bZ   except  -2kF'(.>  )log  r/a(r) 
are  small,  therefore  on  the  envelope 

(38)      log  ^  =  [2kP'(r)]"^  =  0(5**^)    . 

Hence,  at  the  shock,  (10)  may  be  approximated  by  (26)  ar^d  the 
second  term  in  (6)  neglected. 

The  shock  conditions  again  five  that  at  any  point  the  shock 
makes  equal  angles  with  the  characteristics  from  either  side  of 
It;  that  is   U  =  p-(u, +c, +Up+Cp )   where  subscripts  1  and  2  denote 
the  values  on  the  two  sides  of  the  shock.   If  the  shock  is 
t  =  (r-a  )/c   -  gT(r),  and  characteristics  t  =  V^,   V  =■  Z ^ 
{V-y    <  L'2)    intersect  on  it 


-1( 


(UO)      g£(r)  =  kr"^JF(-3^)  +  F{T^)] 

and  from  the  equations  of  the  shock  and  the  characteristics 

(111)  E,^{r)   =  2kF(ri)lo£;  ^  -  r^   , 

=  2kP(r2)lo£  f-  -^2   . 

(rince  the  variation  of  aCr)  Is  small  and  r  Is  relatively  larf,e 
at  the  shock,  a(r)  has  been  set^'  equal  to   a^^  =  aCr  ),  where  T 
is  the  coTiTion  value  of  f,    and  Tp   when  t'ney  are  equal,)   To  solve 
{kC)   and  (Ul)»  the  relation  bet'jeen  ?-,  and  ?'_,  in  order  that 
they  ':neet  on  the  shock,  is  first  fouiido   From  (li-l) 

{k2)      2g|(r)dr  =  2k[?(>^)  +  FCr/g)]  r°^<i'-^ 

+  \2kFHi\)ioi  ^  -  ijdr^^  +  \2kFH^2nor  {-  -  ijdr2 

and  from  (IjC)  the  terms  in  dr  cancel.   From  (Ul),  lo£-  r/a^  may 
be  obtained  in  terms  of  f ,  and  ?,p  as 

m)  i-^—  =•  %-z~-        , 

2k  loc  ~         -2  I 

and  (14.2)  then  becomes 

(Uit)  [F'(r3^)(r2-^i)  -  (f(^2^  ..F(r^))]d?^ 

After  some  manipulation  this  relation  may  be  written  in  a  form 
suitable  for  inter  ration: 

F(r2 )dr2  -  F(>'^)d?^  =  I  d[{r^'.z^)[?{z'^)-^F{r^)]] 

whence 


By  the  argument  following  equation  (30). 
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(U5)    J  P(r)dr  =  |  (r2-ii)[p(r2)+p{?-^)} 


the  constant  of  Integration  belnf  zero  since  '^p  =  'c-,      must  be 
a  solution.  V'lth  the  relation  between  T.    and  f^   known,  r  and  g, 
are  given  by  (lj.3)  and  (i+l)  In  terms  of  ?',  ,  say;  hence  the  shock 
Is  determined. 


Low  these  shock  results  are  exactly  the  same  as  those  en- 
countered and  discussed  in  Pef«  1,  eyceot  that  log  r/a.  appears 

1/2  ' 

in  place  of  r  '  j  the  following  is  taken  practically  word  for 

word  from  T^ef.  1  with  appropriate  changes  in  the  formulae c   The 
results  (14-5)  and  (t4.3)  have  a  very  simple  geometrical  interpreta-- 
tion  on  the  P-curve,   The  left  hand  side  of  (1|5)  is  the  area 
between  the  curve,  the  axis  and  the  two  ordlnates  '£'  =  r^  , 
Z=  t^;   whilst  the  ripht  hand  side  is  the  area  betvreen  the  same 
ordlnates,  the  axis  and  the  straight  segment  joining  the  points 
T    and  Z"p  of  the  curve.   Hence  the  condition  {hS)   requires  that 
the"lobes"  cut  off  on  each  side  of  the  curve  by  this  segtnent 
must  be  equal  in  area.   Moreover  from  (li.3)  the  gradient  of  this 
segment  is  equal  to  (2k  log  r/a.  )"  ,   A  typical  segment  with  the 


Itjbes  shaded  is  shown  in  Figure  1, 

The  path  of  the  shock  may  be  traced  by  considering  the  v.'hole 
series  of  segments;  typical  ones  11*,  22',  33'  are  shovm  in 
Figure  3«   The  limiting  segment  11',  for  which  the  lobe  area 
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has  decreased  to  zero.  Is  the  tangent  to  the  curve  at  the  point 

of  Inflexion  2  =  T'  .   Since  it  has  the  largest  gradient  it  cor- 

o 

responds  to  the  smallest  value  of  r  on  the  shock;  here  the  shock 

starts  with  zero  strength  (since  ?{t)    is  continuous  at  t  =  Z'  ) , 

c 

The  gradients  of  the  successive  sofments  then  decrease,  corres- 
ponding to  increasing  i^,  tending  ultimately  to  be  horizontal  as 
TTp  ->  oo   and  Z\   ->  S'  +  ,  the  zero  cf  F(r)<i   From  (8)  the  pres- 
sure jump  at  the  shock  at  any  point  is 


ike)  p^  [Fir^)   -  PCr^)?  /r 


as  r  increases,  the  shock  first  increases  in  strength  and  then 
decays.   To  study  the  decay  at  lar^  e  distances,  the  approxima-- 
tions  Tp  large,  ?',-?  small  may  be  used.   The  integral  of  F{^) 
from  C  to  oo  exists,  therefore  F{Z)   =  oiZ      )      as  r  ->  oo  ; 
hence  (h5)  i^iay  be  approximated  as 


00 


o 

and  ()'3)  as 

r  .-        2j.^  F(r)dr 

(U8)      ?k  log  ^=  -  -p-21^^ o_         , 

Then  from  (Ul),  g^  =  Z'^i,V^'Z^)F^  {X^)\o^^    r/a,^  -  ^q,  and,  from  (J+8), 
^1  ''^   -^i°o"'^log-^/  r/a^  -  r^,  where 

(U9)     bj  -  -  Xii  /   F(r)dr 

That  is,  the  shock  arrives  at  the  distance  r  at  a  time 
{$0)  b^logl/2  JL 

later  than  the  (asymptotically)  undisturbed  characteristic 

t  =  (r-a  )/c   +  f  .   It  was  found  in  (35)  th-'t,  at  large  dis- 

1/2 
tances,  the  head  shock  was  a  distance   b  lor  '      r/a,  aheqd  of 
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the  same  characteristic.   In  fact  the  shocks  "balance"  since 
b  =  bj^  by  (16). 

It  would  be  expected  that  the  more  jfeneral  procedure 

adopted  for  the  second  shock  could  be  applied  to  the  head  shock, 

and  indeed  this  is  true.   If  the  ^-curve  Is  extended  to  Include 

the  T-axis  for  TT  <  0,  where  the  characteristics  are 

t  =  (r-a  )/c   +T,  'v  =   constant,  segments  cutting,  off  equal 

lobes  of  equal  area  determine  the  shock,  and  the  slope  of  the 

segment  is  (2k   lof(r/a  ))"  .   Such  segTjents  are  shown  as 

aa' ,bb' ,cc' , . . .  in  Figure  3;  the  llmitinr  one  lies  alon,-  the 

discontinuity  at  1^=   0,  and  since  its  rradient  is  Infinite  the 

corresponding  value  of  r  is  a  ,  i.e.  the  shock  is  attached  to 

the  initial  particle  path,   '^he  distance  increases  as  the  sej^ment 

gradients  decrease,  and  as  r  — >  oo  ,  the  segments  tend  to  the 

horizontal  with  TT,  — >  -oo   and  '^p   -^f-*   Thus  Figure  3  gives 

a  simple  interpretation  of  the  flow  pattern.   The  characteristics 

0  <  V  <  V       be  betvjeen  the  shocks:  those  for  which  0  <  ^  <  ^^ 
c  o 

meet  the  head  shock,  and  those  for  which  '5'  <■  X  <  t-      the 
other,  the  asymptotically  straight  characteristic,  Z  =  T  ,    alone 
reaching  infinity.   The  "segments"  give  a  clear  account  of  the 
progress  of  the  shocks  to  infinity. 

At  large  distances  between  the  shocks,  2^  is  near  to  ^'qJ 
therefore,  from  (10), 

(51)    t  =  1:1!°-  -  2kF(r)iog  ^  +  T^      , 

approximately.   Hence,  from  (6), 

/r'0\  r  _       O     _  2y       *  O'    O     O  O 

(52)      i   -  -^ ^  -  :^ --y- 

Y  p^c^   Y     r  log  ^^^ 

At  the  head  shock,  the  value  (from  (35))  is  2( Y+l)"■^br'■^log"■^r/a(-■ 
ln  agreement  i:lth  (3<j)»  and  at  the  second  shock  the  value  is 
minus  this  quantity.   between  the  shoclcs  4  falls  linearly  with 
t  (at  fixed  r)  at  a  rate 

o  c 

^  r  log  -,X_   ' 


.  Oi 


/?.?;..-.*  r,' 


..  .'^'   ::.  ' 
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which  depends  only  on  y*  c  and  the  distance  r;  It  is  inde- 
pendent of  the  particular  explosion  or  Initial  disturbance 
considered. 

i-';ore  complicated  F-curves  tn=in  the  basic  one  in  Figure  3 
may  arise,  and  there  may  be  more  tmn  just  two  shocks.   All 
cases  may  be  described  similar  to  the  above  and  the  details  may 
be  seen  from  the  account  in  Fef.  1.   They  are  not  repeated  here. 
It  is  sufficient  to  remark  that  in  general,  there  are  ultimately 
only  tiijo  shocks  as  r  — >  oo  ,  since  any  others  e'-entu-^lly  run 
into  these  two;  thus  as  r  — >  oo   the  behavior  of  the  flow  is 
eyactly  the  same  as  described  ^bove  for  the  P-curve  of  Figure  3» 
If  PC^)  has  a  discontinuity  at  '^  =  -^  say,  it  mqy  be  shovm  that 
(1'3).  ikS)    (and  hence  the  "segment"  interpretation)  and  (Li6) 
still  apply,  but  near  the  body,  the  expression  (Li-l)  for  g,  is 
modified.   Since  the  "first  segment"  is  vertical  in  this  esse 
the  shock  starts  with  r  =  a,   (as  it  should,  of  course  ,  since 
a^  is  now  the  value  of  a(t. )  at  the  discontinuity  7/  =  P) .  An 
example  with  a  complicated  F-curve,  involving  more  than  two 
shocks,  will  be  given  in  a  later  report,  in  which  the  shocks 
produced  by  the  oscillating  gas  sphere  of  an  underwater  explosion 
will  be  described, 

part  II.   The  Kirkwood-Bethe  Theory  for  the  Propagation 
of  Shocks  in  I'.'ater. 

The  Kirkwood-Bethe  theory  will  be  presented,  as  explained 
in  the  Introduction,  as  a  suggested  method  of  modifying  the  theory 
of  Fart  I  in  order  to  extend  its  application  to  a  non-weak  head 
shock.   The  results  :--hich  vfill  be  obtained  in  this  part  and  the 
arguments  given  below  are  limited  to  propagation  in  t>?ater,  and 
there  is  certainly  no  reason  vjhy  they  should  be  any  improvement 
over  the  results  of  part  I  if  the  medium  is  air.   The  explar^ation 
of  this  is  that  inany  of  the  expressions  and  conditions  used  In 


Since  r  is  lar- e,  the  particular  value  a(T  )  is  not  important; 
any  typical  unit  of  length  In  the  problem  may  be  inserted. 


•r    :'i    jv  :^:. 


19 


Part  I  are  true  not  only  for  weak  shocks  and  disturbances  but 
for  stronc-  ones  provided  that  the  medium  is  water  (this  is  not 
so  for  air).   In  fact  the  method  of  derivation  will  be  to  follow 
through  the  steps  of  Part  I,  making,  only  those  assumptions  which 
seem  reasonable  also  foi  a  non-weak  shock  in  water. 

One  of  the  most  important  advantages  of  v/ater  from  this 
point  of  view  is  that  oven  for  quite  strong,  shocks,  the  change 
of  entropy  is  small."   Kence,  since  entropy  is  conserved  on  the 
particle  paths  behind  the  shock,  the  change  of  entropy  m.ay  be 
neglected.   Therefore  p  may  be  taken  as  a  function  only  of  p, 
the  appropriate  form  being  p  =  K-.p^   «-  K .  ,  with  y  —  7   in  the 
case  of  water.   Then,  introducing  the  velocity  potential  cl>,  the 
momentum  equation  may  be  integrated  to  5ernoulll's  equation 

2  c^ 

(53)      ^r  -^  iu2  +  A  =  °° 


Y-1  •  Z  '^   ■  -t    Y-1 

In  the  theory  of  small  disturbances  of  a  compressible  fluid 
presented  in  Part  I,  r^.    is  constant  on  a  characteristic.   Cn 
the  other  hand,  in  the  eyact  theory  of  an  incompressible  fluid, 
4  =  CT  (t)/r  and  the  characteristics  have  degenerated  into 
lines   t  =  constant,  so  that  again  r^.    is  constant  on  a  charac- 
teristic.  In  view  of  these  two  alternative  cases  the  result 
that   4.  =  -F(r)/r,  where  r  is  the  characteristic  variable,  is 
taken  over  from  part  I,  and  used  ever,  if  the  disturbance  is  not 
weak  and  the  fluid  has  (small)  compressibility.   Then  Bernoulli's 
equation  (53)  becomes 

2   2 
(5U)      ^;^  .  1  u2  =  £fl   . 

To  complete  the  solution  another  expression  involving  u 

and/or  c  is  required.   The  expression  in  Part  I  lor  u(=q  ), 

unlike  that  for  ci.   has  not  the  form  that  u  would  have  in  the 
_5l 

"*"'  Actually  it  is  only  necessary  to  neglect  the  entropy  in  ob- 
taining -Bernoulli's  equation  ((53)  below)  and  in  fact  the 
error  in  (53)  is  still  small  even  when  the  entropy  change 
starts  to  be  appreciable.   In  Ref.  3,  the  error  is  estimated 
to  be  5°/o  when  the  pressure  behind  the  shock  is  20,0CC 
atmospheres . 
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exact  compressible  fluid  theory;  the  first  term  in  (6)  would  no 
longer  appear,  yet  it  was  found  that  this  first  ter.ii  dominates 
the  second  term  in  the  neighborhood  of  the  shock.   Hence  (fc)  is 
not  suitable.   However,  the  previous  remark,  that  only  the  first 
term  in  (6)  was  important  near  the  shock,  leads  to  a  valuable 
suggestion.   This  first  term  in  the  expression  for  u  is 

(55)      cr=  -^  (c-c  ) 

in  the  weak  shock  theory;  we  now  ask  if  u  may  be  approximated  by 
ceven  in  the  neighborhood  of  a  non-weak  shock.   The  shod  con- 
ditions show  that  this  is  so:  values  of  vi/c   and  cr/c   are  shown 

'  '  o        o 

in  the  Appendix  for  values  of  E,   ranging  from  0  to  8o/3 ,  (for 
water  ?  is  more  relevant  than   (p-p  )/p  —  3,000  i^   as  a  measure 
of  the  shock  Intensity),  and  it  may  be  verified  that  the  approxi- 
mation is  good.   !loreover  if  the  problem  were  that  of  plane  waves 
Instead  of  spherical  waves,  it  vjould  be  true  that  u  =  o"  to  the 
same  order  of  approylmation  throughout  the  flow  .   (The  Fiemann 
variable   u  -  2c/y-1,  and  therefore  u-o^,  would  be  constant  on 
each  "negative"  characteristic;  since  u-cr  is  approximately  ?;ero 
at  the  shock,  it  would  follow  that  \x-6'  is  zero  to  this  order  of 
approximation  everywhere,)   Thus,  even  for  non-weak  shocks, 
u  —  o~  at  the  shock  and  also  for  some  distance  behind  (since  the 
spherical  nature  of  the  flow  is  not  apparent  immediately  and  it 
behaves,  at  first,  like  one -dimensional  flow).   But,  one  important 
point,  which  is  easily  overlooked,  remains:   In  determining  the 
characteristic  variable  2- (which  is  used  in  finding  the  shock), 
it  is  necessary  to  integrate  the  equation  dt/dr  =  (u+c)"   along 
the  whole  charactei istic ,  even  in  regions  away  from  the  shock 
where  it  is  not  true  that  u  —  O'.      Kow,  in  Part  I,  the  same  thing 
occurred,  but  it  v;as  still  found  that  the  contribution  to  this 
integration  from  the  second  term  in  (6)  (the  deviation  of  u  from 
c  in  that  case)  could  be  neglected  in  determining  the  shock, 
although  it  would  be  important  in  other  regions.   (This  was  the 
result  deduced  after  (214.).)   It  is  assumed  that  this  is  still 
true  in  the  case  of  a  non-weak  shock.   The  author  feels  that  this 
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assumption,  that  the  Integrated  effect  of  u-cT  is  nef ligible  at 
the  shock,  is  the  only  one  in  the  found q.tion  of  the  Klrkwood- 
Bethe  theory  which  can  be  seriously  questioned;  it  is  the  only 
one  taken  over  from  part  I  which  has  no  additional  support  from 
"non-weak  results."  Igevertheless  its  use,  \intil  it  is  refuted 
by  experiment  or  numerxcal  calculations  or  supplanted  by  a  better 
theory,  seems  very  reasonable,, 

The  solution  for  cr(r,r),  obtained  from  (5U)  and  (55)  with 
u  =  cf "  is 

(56)  c^o'd  +  Pd-)  -  I^   , 
where 

(57)  p  =  kc^  =  ^  (y+D/Cq   . 

The  characteristic  variable  t-'(r,t)  occurring  here,  is  found  from 
the  condition  that  on  f  =  constant,  dr/dt  =  u+c.   Ivith  u  =  o', 
(55)  gives   u+c  =  c  (l+2po);  therefore 

(58)  t  =    r         _T_^^r^  ^  ^       ^ 


^(^j   c^'il+2P(3-(r,r)j 


o 


being  defined  uniquely  as  before  to  be  the  value  of  t  when 
r  =  a.   Equations  (56)  and  (57),  giving  O"  and  t  as  functions  of 
r  and  the  para.Tieter  '^ ,   provide  the  solution  to  the  flow;  u  is 
approximately  equal  to  d  near  the  shock,  and  p  is  given  In  terms 
of  crby 

o  o   ^  o-  \         o/ 


p. 


The  arbiti'ary  function  F(T)  depends  on  the  boundary  conditions 
on  r  =  a(T).   In  this  case  F{t)    cannot  be  determined  from  equat- 
ing u  to  a' (T)  on   r  =  a(f),  since  an  accurate  expression  for  u 
in  this  rep-ion  is  not  known;  but  ifj  for  example,  p  is  given  as  a 
function  of  r,  P("')  can  be  found  from  (59)  and  (56)  with 
r  =  a  ( ^' ) . 
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The  shock  relations,  (22)  and  (23),  which  were  u?ed  for 
weak  shocks,  are  ^^ood  approxlnatlons  also  for  non-weak  shocks  In 
water  (see  Appendix);  (22)  Is  already  satisfied  by  taking  u  =  o^ 
and  the  other,  which  may  be  written 

(60)      U  =  Cq(1  -^  Po-)    , 

is  now  used  to  determine  the  shock.   First,  it  is  convenient  to 
express  iS'o)   differently.   ProTi  (56), 

therefore 

Where  c(T)  is  a  known  function  of  T ,  being  the  value  of  o"  on 
r  =  a(cr)  i.e.  the  solution  of  (56)  with  r  replaced  by  a(7r). 
Hence  (5&)  becomes 


(61)      t  =  r  -  ^^  |A(a)  -  B(r)] 


% 


where 

r 


^  ^     ^  ^ '  J  ;g(i.p<,)^   ^^  ^°<~^j "  oTiwj-   ' 

3(r)  =  A(c(r)) 

In  Part  I,  it  was  convenient  to  determine  the  f auction  r  =  B{V) , 
giving  the  distance  r  of  the  shock  on  each  characteristic;  in 
the  present  case  it  is  more  convenient  to  find  the  value  of  cf  at 
the  shock  as  a  function  of?".   Let  this  function  cr(R(r),T)  oe 
denoted  by  8(2' )♦  'hen  b{T)    is  determined,  r  is  given  from  (56) 
as  R(r),  where 


(63)      R(r)  = 


F(^) 


c^s(r)>i+ps(r)? 

o     (.        ^ 


:v:  ■■& 


;  J-«' 


T::--- 


1-  , 
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and  t  Is  given  from  (61)  as  a  function  T(?),  where 


(61+) 


Tcr)  =  r-  ^^^UC3{l:))  -  B(r)] 


; 


Thus  the  shock  and  the  values  of  cr,    u,  p  at  the  shock  will  be 
known  with  T  as  parameter. 

The  shock  velocity  Is  R»(T)/T'(T);  therefore.  In  ( 60 ) , 


(65) 


P'(n  =  c^(l+P3)T'(r) 


substituting  from  (63)  and  (6I4.),  (65)  gives  the  following  first 
order  equation  for  3(f) : 


(66)  --^ 


ppcn 


c^s(l+p3) 


d-3        2(3   r.         l+ps  3+2ps      ,   r^,/.-^N 


-  1    +     1      ^ 

-  1   +  -^ 

^o 


2(l+ps)-| 

^^-^F(r)3(r)] 


As  for  the  corresponding  equation  (26)  of  Fart  I,  (66)  may  be 
obtained  in  a  form  suitabls  for  integration  by  multiplying  by 
P(T);  on  integration 


(67) 


P 
°o 


.^^rA^os^   -^^J 


C   '  "^^    2(l+psi? 


.t 


-  o 


=  I  F(rMdr«  + -^  / 'F('c')-^(?(r')3(t'))dr' 

-'o  C   .y'  o 

o 

This  formula  determines  3(r)  and,  therefore,  the  shock.   It  is 
convenient,  however,  to  simplify  the  second  Integral  on  the 
right  hand  side.   After  some  manipulation,  it  is  found  that 

.2- 

(68) 

c 
o 


I  =  \      F(r«)  g^T  (F(r')B(r'))cir' 

c  <J  o 


=  -^SB(T)P^(r)-B(0)F^(0)^  +  -^  i  F^(r')  ^|£^  dr 
2c„  C  J    2c^  -^o 


2U 


then,  since  dS/dT  =  a~^(  l+!3a)"^da/d:%  from  (62),  and 

(69)  c^a(l+pa)  =  P(V;)/a(r) 
I  Is  found  to  be 

(70)  I  =  -^   ^B(r)P^(r)-B(0)P^(C)?  ■"  I  /    '  a2(ir')da(r') 

where   a  =  a(C).   Using  this  expression,  with  the  value  for 
BC?")  substituted  and  the  last  term  Integrated  by  parts,  and 
remembering  that   3(0)  =  a   since  at  T  =  0   the  shock  meets 
the  particle  path,  (67)  becomes 

( 


-,n  N      3F^(r)    ,^(-,         l-iPs  3+2(33  7         S,         1+Pa  3+2pa   7   ) 

71)     Eil-^^      l^iog  -p-L-  -  ^i^^5    -  [log  -^  -  2TrT(J5T5/ 


o- 


=  /  P(t')dr  ~  /  a(r')a(r')  ^^fTP  dn   . 

Jo  ^o  '' 

Klrkwood  and  Bethe  arrive  at  this  result  (although  not 
written  In  the  same  form)  but  then  remark  that  for  most  purposes 
certain  terms  may  be  omitted.   In  the  notation  and  derivation 
adopted  here,  their  suggestion  amounts  to  approximating  I  (equa- 
tion (70))  by 

c. 


(72)      J  =  1  ^1|I  j^p2(->;)  _  p2(o)^    ; 


it  is  difficult  to  believe  that  this  is  a  valid  approximation. 
For,  when  T  is  small,  J-^  c^^B( 0)P( 0)P« (0 )r,  whereas 

(73)      I~-  j  c^^B(0)P(0)P'(0)+  I  c;^F^(0)B'(0)+  ^  a^a'(0)|r 

=  S  c^^B(0)P(0)P'(0)+a^a'(0)  j'?' 

from  (62)  and  (69);  when  P(r)  ->  0   (corresponding  to  large 

1-2     2 
distances  on  the  shock)   J  ->  -  2*  °o  ^^^^^  (*^^*  whereas 
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ilk)  I  ->  -  I  c;^3(C)F^(C)  +  J  /      a^(?:Mda(J') 


o 


Now  In  the  special  case  of  a  small  disturbance  with  a  7/  0   the 
dominant  term  In  BCT)  Is  -l/a  and  0   a{T)    =  F{T)/a{Z')      so  that 


-2„,..^,..„....  _   .2 


o 


c^'^B(0)P(C)P' (0)  =  -a^a«(0)(l+C(e))  ,  where   e  =  a'(0)/c^;  hence 

the  additional  term  In  (73)  actually  cancels  the  term  retained 

1-2    -2      12 
In  J.   Again,  for  this  case,  -  2"  °o  3(C)F^(C)  =  ^   ^0%'  therefore 

the  additional  term  in  (7I4-)  is  certainly  not  negligible.   Thus 

it  appears  that  I  should  not  be  approximated  by  J. 

The  shock  relation  (71)  has  been  written  in  a  form  which  is 
especially  convenient  for  comparison  with  the  corresponding 
relation  (30)  of  Fart  I.   It  is  seen  immediately  that  in  the  case 
of  weak  disturbances  with  s  and  a  small,  and   c  a  =  P(Z")/a, 
c  3  =  F(r)/P,  the  left  hand  side  of  (71)  becomes 

—2  -2 

Pc   log(a/s)  =  pc~  log  R/a;  also  the  second  term  on  the  right 

hand  side  becomes  (/■    PCl' )c"*  a'(^')d'-',  and  since  a'(f)/c   is 

J  o  ^sc 

small  this  term  may  be  neglected  in  comparison  with  '<      F(^')d'^'. 

*-  o 
Hence  (30)  is  obtained. 

Finally,  consider  the  propagation  at  large  distances. 
Since   s  -c>  0,  it  is  clear  from  (71)  that  Z   ~>  X    ,   where  7r^   Is 
the  first  zero  of  P(^);  hence  3(T)  and  a(r)  are  small  and  (71) 
becomes 

?  ?  f  "^  o  r  ^o 

(75)   Pc"'^p''(r)log  R/a  =  j    F(r')dr'  -  ;    a(r')a(2'')a'(r')cir<   . 

i>  o  *>  o 


The  only  change  to  the  asymptotic  formulae  of  part  I,  therefore, 

is  in  the  value  of  the  constant  b,  In  which   /  "  F{Z*)dZ^      must 

^o 
be  replaced  by  the  right  hand  side  of  (75)«   The  variation  in 

pressure  behind  the  shock  is  apain  given  by  iSZ^i . 

Thus  the  extension  of  the  results  of  Part  I  is  quite  simple; 
how  effective  it  is  must  be  decided  by  numerical  work  on  special 
cases. 
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Appendix 

Nvunerical  values  of  u,  C,    U   for  a  shock  in  water  are  given, 
for  various  values  of  the  pressure  Jump,  in  Ref.  3«   Some  of  the 
results  are  set  out  here  in  order  that  the  relations  u  —  0"  and 
U  —  c  (1+Pcr)   at  the  shock,  whichwere  used  in  the  Kirkwood- 
Bethe  theory  may  be  verified.   The  values  taken  for  the  constants 
Y  and  c   are  7.l5  and  II4.65  metres/sec,  respectively;  K  varies 
slowly  with  the  entropy,  its  value  ahead  of  the  shock  being 
3.0i4-7  kilobars.   The  pressure  is  measured  in  kilobars  and  it 
should  be  remembered  that  1  kllobar  is  approximately  ICOO  atmos- 
pheres, so  that  the  strength  (p-p  )/p   is  approximately  1,000 
times  the  numbers  appearing  in  the  pressure  column  and  4  is 
approximately  1/3  of  these  numbers. 


p 

u 

cr 

u 

c^d+po-) 

kilobars 

m/sec 

m/sec 

m/sec 

m/sec 

10 

i;26 

I4.I6 

2335 

23II1 

20 

669 

663 

2880 

2817 

30 

898 

659 

3320 

3217 

ko 

1075 

1030 

3690 

3566 

50 

1235 

1165 

i+020 

3682 

60 

1380 

1330 

1+325 

1+176 

70 

1510 

11+65 

1+6 10 

1+1+51^. 

80 

1625 

1600 

1+885 

Ii729 
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